Nuclear magnetic octupole moment and the hyperfine structure 

of the 5/^3/2,5/2 states of the Ba + ion 

K. Beloy, 1 A. Derevianko, 1 ' 2 V. A. Dzuba, 2 G. T. Howell, 3 B. B. Blinov, 3 and E. N. Fortson 3 

1 Physics Department, University of Nevada, Reno, Nevada 89557 
2 School of Physics, University of New South Wales, Sydney 2052, Australia 
3 Department of Physics, University of Washington, Seattle, Washington 98195 

(Dated: April 28, 2008) 

Abstract 

The hyperfine structure of the long-lived 5D^/ 2 an d 5Z?5/2 levels of Ba + ion is analyzed. A 
procedure for extracting relatively unexplored nuclear magnetic moments Q is presented. The 
relevant electronic matrix elements are computed in the framework of the ab initio relativistic 
many-body perturbation theory. Both the first- and the second-order (in the hyperfine interaction) 
corrections to the energy levels are analyzed. It is shown that a simultaneous measurement of the 
hyperfine structure of the entire 5Dj fine-structure manifold allows one to extract without 
contamination from the second-order corrections. Measurements to the required accuracy should 
be possible with a single trapped barium ion using sensitive techniques already demonstrated in 
Ba + experiments. 

PACS numbers: 32.10.Fn,31.15.am,21.10.Ky,27.60.+j 
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I. INTRODUCTION 



A nucleus as a source of the electro-magnetic fields is conventionally described using a 
hierarchy of static electromagnetic moments: magnetic dipole (Ml,/x), electric quadrupole 
(E2, Q), magnetic octupole (M3,f2), etc. Interaction of atomic electrons with these moments 
leads to the hyperfine structure (HFS) of the atomic energy levels. While the first two 
moments, \x and Q, have been studied extensively, the octupole moments remain relatively 
unexplored. 

While octupole moments may be approximated using the nuclear-shell model [lj, the 
correct values depend strongly on nuclear many-body effects and, in particular, on core- 
polarization mediated by the nucleon spin-spin interaction |2J . Senkov and Dmitriev [2] car- 
ried out a nuclear-structure calculation of Q for 209 Bi. In this particular case, the polarization 

n 

effects enhance the shell-model values by a factor of three. According to Dmitriev [3J, a sys- 
tematic study of octupole moments will help place constraints on the poorly-known isoscalar 
part of nuclear spin-spin forces. In another, even more striking example, the deduced value 
of f2 of 133 Cs is forty times larger than the shell-model value [4| (this has not been analyzed 
yet in nuclear theory). We note that in the case of Cs, the shell-model value is strongly 
suppressed due to an accidental cancellation between the orbital and spin contributions of 
the valence proton to the magnetic-octupole moment. 

Measuring the effects of the octupole moments on the hyperfine structure was so far lim- 
ited to a small number of atoms: CI fl, Ga 0, Br fl, In fl, V fl, Eu Q, Lu 11], Hf 12], 
13] . Since deducing Q from a measurement requires knowing atomic-structure cou- 



and Bi 



plings, previous analysis focused primarily on isotopic ratios because the electronic coupling 
factor cancels out when ratios of HFS constants are formed. An exception is the measure- 

ri 

ment on the 6P3/2 state of 133 Cs |4|, where sufficiently accurate calculations are possible. 
In a recent paper [14j, we argued that an accurate deduction of the octupole moments is 
feasible for metastable 3 P2 states of alkaline-earth atoms. 

Ba + , being an atomic system with one valence electron outside a closed-shell core, also 
presents a case where both high-accuracy measurements and high-accuracy calculations are 
possible. The goal of this paper is to analyze the hyperfine structure of the 5D 3 / 2 and 
5-D5/2 levels of Ba + , and to show that Ba + is a particularly favorable case for measuring 
octupole moments, for both theoretical and experimental reasons. Both 5D levels belong 
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to the same fine-structure manifold. We take advantage of a simultaneous analysis of the 
hyperfine structure of both levels and show that such an analysis allows to eliminate the 
potentially troublesome second-order hyperfine electron structure term thus providing a 
powerful consistency test for the measurements and calculations and possibly improving the 
accuracy of deducting Q. Furthermore, 5D 3 / 2 and 5D 5 / 2 are each long lived metastable 
states in Ba + , with lifetimes of about 80s and 30s respectively. Hyperfine intervals thus 
could be measured in principle to well below 0.1Hz, much better than in previous octupole 
experiments; techniques to exploit the sensitivity inherent in the 5D levels have already been 



developed using single trapped Ba + [15 ]. 



Barium has two stable odd isotopes, 135 and 137. Both isotopes possess nuclear ground 
states with spin of / = 3/2 and of positive parity. The values of the magnetic and quadrupole 
moments are 

135 Ba :/i = 0.837943 li N , Q = +0.160 b, 

137 Ba : li = 0.937365 li n ,Q = +0.245 b , (1) 

where lin is the nuclear magneton and b (barn) = 10~ 24 cm 2 . Both isotopes have unpaired 
neutrons in the single-particle state, and from the single-particle (shell) model [l| we 
may estimate the octupole moment to be 

tt sp = 0.164/i W (r 2 ) « 0.0385 ll n x b , (2) 

where we used the rms value of the nuclear radius (r 2 ) 1 / 2 = 4.84 fm. Of course, the shell 
model is only an approximation. For example, in contrast to the known properties ([1]) for 
the two isotopes, this model would produce identical values of the magnetic moments and 
vanishing values of the quadrupole moments. 

This paper is organized as follows. First we recapitulate the theory of the hyperfine 
structure of atomic levels, including octupole moments and the second-order effects. We 
present specific formulae for the Ba case and show how to extract the octupole constant 
C from hyperfine intervals (to be measured). Further we compute the electronic structure 
factor required for extracting the octupole moment from the constant C . Finally, we present 
a brief description of an experimental method with single Ba + to determine the octupole 
moments of 135 Ba and 137 Ba. We follow the notation and formalism of the recent paper 14]. 
Unless specified otherwise, atomic units, % = |e| = m e = 1, and Gaussian electromagnetic 
units are employed throughout. 



II. PROBLEM SETUP 



Hyperfine-interaction (HFI) Hamiltonian, describing coupling of electrons to various nu- 
clear moments may be represented in the tensorial form of 

#HFI = ^2 ' T k,q T k,-q- 

k,q 

Here the spherical tensors (of rank k) T? act on the electronic coordinates. Tensor op- 
erators Tj? q are the components of the nuclear electric and magnetic 2 fc -pole (MJ and EJ; 
T,P-even) moment operators. In particular, the conventionally-defined magnetic-dipole, 
electric-quadrupole, and magnetic-octupole moments of the nucleus are proportional to the 
expectation values of the zero- component (q = 0) operators in the nuclear stretched states 
|J,M/ = I): fj, = (T fc n =1 )/, Q = 2(T fc n =2 ) 7 , and Q = -(T fc n =3 )/. Explicit expressions for 
the electronic operators and the corresponding reduced matrix elements are tabulated in 
Ref. Q. 

The conserved angular momentum F for the hyperfine coupling is composed from atomic, 
J, and nuclear, I, angular momenta: F = I + J. It is convenient to work in a basis spanned 
by the eigenfunctions \ jIJFMp) which is formed by coupling atomic, \jJMj), and nuclear, 
\IMj), wave functions. Here 7 encapsulates remaining electronic quantum numbers. For 
1 = 3/2 each of the 5Dj levels splits into four hyperfine components: 5D 3 / 2 has hyperfine 
components F = 0, 1, 2, 3, and 5D 5 / 2 has components F = 1,2, 3, 4. 

Owing to the HFI's rotational invariance, a matrix element of the HFI in the \ jIJFMf) 
basis is diagonal in the quantum numbers F and Mp. If we limit our system of levels to 
only the 5Dj fine-structure manifold, the hyperfine components F = 1,2,3 of the 5D 3 / 2 
and 5D 5 / 2 levels become coupled. The intervals within each manifold may be parameterized 
using the conventional hyperfine constants A, B, C and the second-order corrections (in 
HFI) i] and (. Constants A, B, C are proportional to nuclear moments fi, Q, and Q. Ml- 
Ml correction rj is of the second order in /i, and £ comes from a cross-term between Ml and 
E2 parts of the HFI. Second-order corrections are suppressed by a large energy denominator 
equal to the fine-structure splitting between the 5Dj levels. 

The energy intervals SW F J ^ = W F JS> — W F J \ within each fine-structure manifold 5Dj are 
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as follows. For 5_D 3 / 2 , 



8W[ m = -2A + B + 28C + ^r], 
1 75 

,(3/2) „ , „ „„ 1 1/3 



W -' = -3A- B -8C + -„ + -^C, (3) 



and for 5D 5 / 2 , 



^ = ^-! fl -f c+ ^ + ^- (4) 

In the above equations the HFS constants A, B, and C are all specific to the state of 
consideration while r] and ( represent the same second-order HFS constant. 

If we assume all other second- and higher-order effects are negligible (see justification in 
Section HVl) . then we may solve for the HFS constants A, B, and C in terms of the HFS 
intervals and these two second-order constants. Specifically, solving for the C constants: 

V 31 ' 80 100 1 400 2 2000^ 

C(5D 5/2 ) = --5W[ b/2) + -5W 2 {5/2) - — <W 3 (5/2) + -^=C- (5) 

V 5/ ; 40 1 35 2 112 3 200^ 

The C constants do not depend on the Ml-Ml 77 correction, as was proven in Ref. [ijj] on 
general grounds. 

It is possible to use Eqs. (j^D to cancel the constant £ and therefore eliminate the second- 
order effects from the problem altogether. In doing so, we obtain the equation 

C(5D 3/2 ) + ±C(5D 5/2 ) = -^ (3/2) + T^^i 3 ^ - WV W2 ' 2) 

1 —(5/2) , J_ r M/ (5/2) 1 r W (5/2) 



400 1 350 2 1120 3 v ; 

Since each of the constants C is proportional to the same octupole moment, knowing the 
hyperfine splitting inside each of the fine-structure manifolds provides direct access to Q. 
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III. ELECTRONIC-STRUCTURE FACTORS 



Provided that the measurements of the HFS intervals are carried out, one could extract 
the octupole moment by computing the matrix elements of the electronic coupling tensor 
T 3 e . Specifically, 

/ J 3 j\ 

C(jJ) = -Q 7 J IJ 7J). 

I -J J J 

We may also compute the second-order M1-E2 HFS correction ( by computing the off- 
diagonal matrix elements of the electronic coupling tensors Tf and T%. For / = 3/2, ( is 
given by 

_ 20 ^Q(5D 3/2 \\T^\\5D 5/2 )(5D 5/2 \\T^\\5D 3/2 ) 

\/3 E 5D 5/2 ~ E, 



11 -'5D3/2 

Matrix elements of the electronic tensors are given in Ref. 1141. 

To calculate the electronic-structure factors for the HFS constants C(5D 3 / 2 ) and C(5D 5 / 2 ) 
and for the second-order term ( we employ the correlation potentia. 



order correlation correction operator S(°°) as suggested in Refs. 



17 



method [r| using all- 



181 ] . The method was 



used for Ba + previously [19j for accurate calculation of the parity non-conservation. It is 
also known that the method produces accurate results for the magnetic dipole hyperfine 
structure constants of alkali atoms (see, e.g. Ref. [201]). 

Calculations start from the Hartree-Fock procedure for the closed-shell Ba 2+ ion 

(H - e c )^ c = 0, (7) 

where H is the single-electron relativistic Hartree-Fock (RHF) Hamiltonian 

H = ca-p + (/3- l)mc 2 - Ze 2 /r + V core , (8) 

index c in Eq.Q numerates core states, and V core is the sum of the direct and exchange 
self-consistent potential created by Z — 2 core electrons. 

States of the external electron are calculated using the equation 

(H + t - e v )i) v = 0, (9) 

which differ from the equation for the core (J7J) by an extra operator E. The so-called 
correlation operator E is defined in such a way that in the lowest order the correlation 
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correction to the energy of the external electron is given as an expectation value of the £ 
operator 

8e v = <V«|S[V«>- (10) 

The correlation potential £ is a non-local operator which is treated in the Hartree-Fock-like 
equations (jSJ) the same way as a non-local exchange potential. Solving these equations we 
get the energies and the orbitals which include correlations. These orbitals are usually called 
Brueckner orbitals. 

We use the Feynman diagram technique 17J and B- spline basis set 21] to calculate £. 



The many-body perturbation theory (MBPT) expansion for E starts from the second-order 
and has the corresponding notation £( 2 ). However, we also include two dominating classes of 
the higher-order diagrams into the calculation of E, as described in Ref. 17]. These higher- 
order effects are: (1) screening of Coulomb interaction between core and valence electrons by 
other core electrons and (2) an interaction between an electron excited from atomic core and 
the hole in the core created by this excitation. Both these effects are included in all orders 
and corresponding £ is called Another class of higher-order correlations is included 

in all orders when the equations (jSJ) are iterated for the valence states. These higher-order 
effects are proportional to (E) 2 , (S) 3 , etc. 

The effect of the second- and higher-order correlations on the energies of Ba + is illustrated 
by the data in Table [H As one can see, the inclusion of the correlations leads to systematic 
improvement in the accuracy for the energies. Note that since solving the equations (jUJ) 
for Brueckner orbitals produces not only the energies but also the wave functions of the 
external electron, the better accuracy for the energy should translate into better accuracy 
for the wave function and for the matrix elements. Therefore, we can try to improve the 
wave function even further by fitting the energies to the experimental values by rescaling 
the E operator in Eq. (jUJ). This is done by replacing E by /E, where rescaling parameter / 
is chosen to fit experimental energies. The values of / for different states of Ba + are listed 
in Table [H 

To calculate the HFS constants we need to include extra fields which are the fields of 
the nuclear P-even electromagnetic moments such as magnetic dipole, electric quadrupole, 
etc. This is done in the self-consistent way similar to the RHF calculations for the energies 
in the frameworks of the well-known random-phase approximation (RPA). Corresponding 
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TABLE I: Removal energies of Ba + in different approximations [cm x ] 



State 


J RHF YS l > T, (co > fitted a 


Expt. 6 


6S 


1/2 75339 82227 80812 80685 


80687 


6P 


1/2 57265 61129 60584 60442 


60425 


6P 


3/2 55873 59351 58863 58735 


58734 


5D 


3/2 68138 77123 76380 75816 


75813 


5D 


5/2 67664 76186 75543 75004 


75012 



a /(6s)=0.97^ 
h NIST, Ref. 



22] 



(6p)=0.960, /(5d)=0.934 



equations have the form 

(H - e c )# c = —(F + 5V core )^ c , (11) 

where F is the operator of external field, 5ip c is the correction to the core state due to 
the effect of external field and SV core is the correction to the self-consistent Hartree-Fock 
potential due to the change in field-perturbed core states. The RPA equations ( fill can be 
considered as a linearized (in external field) expansion of the RHF equations ([7]); these are 
also solved self-consistently for all the core states. This corresponds to the inclusion of the 
so-called core polarization (CP) effect. Matrix elements for states of the external electron 
are given by 

(ipvlF+SVcorellPv). (12) 

Dominant correlations are included by simply using the Brueckner orbitals as the wave 
functions ip v in (fl~2]) . There are, however, correlation corrections to the matrix elements 



which are not included into ffl~2l . These are the structure radiation (SR) and the effect of 

n 

normalization of the many-electron wave function [16| . Structure radiation can be described 
as a contribution due to the change in S caused by the effect of the external field: 

(ip v \St\ipv). (13) 

We calculate SR and renormalization contributions using the MBPT similar to the third- 
order calculations presented in Ref. [23|] (second-order in Coulomb interaction and first-order 
in external field). However, we use the "dressed" operators of the external field: F + 5V C 
rather than just F as in Ref. 



231 ] . Therefore, core polarization effect is included in all orders 
8 



in the SR and renormalization calculations. We also use two different basis sets of single- 
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25J, and another is the 



electron states. One is the dual kinetic-balance basis (DKB) set 
B-spline basis set developed at the University of Notre Dame [21 ]. 

The results of the calculations are presented in Table [TT1 Here the RHF approximation 
corresponds to the (ip^- F \F\ip^ F ) matrix elements with the Hartree-Fock orbitals ip^ F . RPA 
approximation corresponds to the (ip^ F \F + 5V core \ip^ F ) matrix elements. Brueckner + CP 
approximation corresponds to the (ip^ T \F + 5V core \tp® r ) matrix elements with Brueckner 
orbitals ip® r , etc. The values of the SR and renormalization corrections listed in Table HT1 
were computed with the DKB basis set. 

TABLE II: Magnetic octupole hyperfine structure constant C of the 5-D 3 / 2 and 5_D 5 / 2 states of 
Ba + and off-diagonal matrix elements of the magnetic dipole and electric quadrupole operators in 
different approximations. 



Approximation 


C(5Ay 2 ) 


C(5D 5/2 ) 


(5D 3/2 \\Ti e) \\5D 5/2 ) 


(5£> 5/2 ||T 2 (e) | |5Z> 3/2 ) 




kHz/(n/(/x N x 


b)) kHz/(n/(/x N x b)) 


MHz //i 7v 


MHz/b 


RHF 


-0.4294 


-0.1514 


-95 


180 


RPA = RHF + CP 


-0.5843 


0.9636 


-1360 


184 


£( 2 > + CP 


-0.6863 


0.9254 


-1496 


222 


S(oo) + C p 


-0.6822 


0.9244 


-1489 


220 


Energy fitting (Br) 


-0.6758 


0.9282 


-1481 


218 


SR 


0.0842 


-0.8472 


280 


14 


Norm 


0.0178 


-0.0287 


42 


-5 


Total (Br+SR+Norm) 


-0.5738 


0.0523 


-1160 


227 



IV. EXTRACTING OCTUPOLE MOMENT 

Our final results for the magnetic octupole hyperfine structure constants are 

^(51)3/2) = -0.585(11) (j^—^j kHz ; 

C(5A, /2 ) = 0.036(16) (— 5— J kHz. (14) 

\fj, N x by 
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Here central values and the errors are found from the scattering of the results due to effects 
of energy fitting and change of basis for the SR and renormalization calculations. Notice that 
the error bars are purely theoretical and reflect the fact that only certain classes of diagrams 
are included in the calculations. In particular, there are strong cancelations between various 
contributions to the C(5D 5 / 2 ) constant, leading to a large, 45%, uncertainty in the value of 
this constant. 

The above error estimates are consistent with the general trend for the experimentally 
known constants A and B of the 5Dj states [26f . Our employed method is off by as much 
as 10% for A(5D 3 / 2 ) and 30% for A(5D 5 /2). The computed values of B generally agree at 
the level of a few per cent with the experiment. Another insight comes from understanding 
that the theoretical method includes the RPA and Brueckner chains to "all orders"; as long 
as these classes of diagrams dominate, the theoretical accuracy is excellent. By contrast, 
in case of the C(5D 5 / 2 ) the result is accumulated due to remaining SR and Norm diagrams 
(see Table HI]), which are computed nominally in the third-order MBPT only. This explains 
the relatively poor accuracy for the 5D§/2 states. As demonstrated in Ref. 27] (at least for 
the constants A and B) the theoretical accuracy could be improved to 1% by employing the 
relativistic coupled-cluster method. 

The results, Eqs. (THl) . may be used to extract the values of the nuclear magnetic octupole 
moment from the measurements. For example, if equation ([6]) is used then 

C(5£>3/ 2 ) + ^C(5D 5/2 ) = -0.581(13) ( — 5— ) kHz. (15) 
10 \Vn x by 

Alternatively, one can use the first equation of (jSJ). Then the correction due to the second- 
order term ( needs to be taken into account. With the values of the magnetic dipole and 
electric quadrupole HFS matrix elements presented in the last columns of Table [III this 
correction reads 

1 f 1.84 Hz 135 Ba+ 

AC(5D 3/2 ) = =( = { (16) 

2000V/2T 3.17 Hz 137 Ba + 



Notice that in this second scenario we advocate using the 5D 3 / 2 hyperfine manifold for ex- 
tracting the nuclear octupole moment because of the poor theoretical accuracy of computing 
electronic couplings for the 5D 5 / 2 state. 

We may evaluate the relative influence of Q on the HFS by using the single-particle 
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(shell-model) estimate for the nuclear octupole moment, Eq. (J2J); we arrive at 



C(5 J D 3/2 ) s - p - « -23 Hz, 
C(5D 5/2 y*- « 1.4 Hz. 

We see that the second-order correction, AC(5-D 3 / 2 ), is below the anticipated value of the 
constant. 

At this point we briefly discuss the effect of all other second- and higher-order terms 
past rj and (, which until this point have been assumed negligible. One might expect 
other second-order dipole-dipole terms which mix in states outside of the 5D fine structure 
manifold to have an appreciable effect on the hyperfine structure. However, the proof in 
Ref. [bj can easily be generalized to show that no second-order dipole-dipole terms enter 
into the equations for the C constants, Eqs. (jSJ). Furthermore, it is found that the leading 
third-order term, the dip ole- dipole-dipole term mixing the fine structure levels, drops out of 
Eq. (jSJ) along with (. Therefore, we can expect the largest terms neglected from Eq. to 
be the second-order dipole-octupole and quadrupole-quadrupole terms; we have estimated 
these effects to both be at the ~ 10~ 3 Hz level. This provides sufficient confirmation that 
all second- and higher-order terms may be neglected in our proposed scheme of extracting 
the octupole moment. 

Finally, using the single-particle approximation for the nuclear octupole moment, we 
obtain an estimated value for the left hand side of Eq. (jSJ) of —22 Hz. Assuming a conservative 
value of 1 Hz uncertainty in the 5D HFS intervals yields an overall uncertainty of 0.017 Hz in 
the right hand side of Eq. (jHJ). In the next section we describe how such measurements should 
be capable of much smaller uncertainties, 0.1 Hz or better. Thus we may conclude that HFS 
interval measurements with readily attainable accuracy, combined with the theoretical result 
( I15p . would be capable of extracting an octupole moment of the estimated size for the 135 Ba 
and 137 Ba nuclei. 

V. EXPERIMENTAL POSSIBILITIES 

The measurements can be carried out by a technique similar to one already used to study 
transitions among sublevels of the 5D 3 / 2 state of Ba + [ljj], in which optical pumping is 
used to place the ion in a particular sublevel, and an RF transition to another sublevel is 
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FIG. 1: The lowest S, P, and D states of Ba , showing the cooling 493 nm and cleanup 650 nm 
transitions plus the 2051 nm and 1762 nm E2 transitions to the metastable D states. 



(a) — ^ 

1762nm / 


(b> RF tJ- F = 1 


(c) -2- ^ 

1762nm / 


(d) 6P1/2-7-F 

I \.650nm 

493nm/ 


F 1 6S " 2 






6S„ 2 



FIG. 2: Procedure for measuring the hyperfine intervals in the 5D^/ 2 state, (a) First, a 1762 nm 
resonant laser pulse transfers the ion into the F = 2 hyperfine sublevel of the 5D5/2 state, (b) The 
RF is then applied to drive a hyperfine transition in the 5-D 5 / 2 state, (c) The second pulse of the 
1762 nm laser depopulates the 5D§/ 2 F = 2 level, (d) To determine if the hyperfine transition in 
the 5-D 5 / 2 state occurred, the cooling and the cleanup lasers are turned on and any ion fluorescence 
is detected. Absence of fluorescence indicates that the RF is on resonance and caused the transition 
to take place. 

detected by the effect of "shelving" as described below. Such measurements are performed 
on a single ion held by radio frequency electric fields in a 3-dimensional effective potential 
well typically ~ lOOeV deep, with the ion at the bottom of the well after being laser cooled 
to a temperature = 10~ 3 degK, with an orbital diameter = 10 -2 /zm. 

The electronic energies of the lowest S, P, and D states of Ba + are shown in Fig. [lj 
The cooling laser operates on the QS1/2 — 6P1/2 allowed El absorption line near 493 nm, 
mistuned slightly to the red of resonance to effect Doppler cooling. A 'cleanup' laser beam 
operates at the 6P1/2 — 5D%/2 transition near 650 nm to keep the ion from getting stuck in 
the metastable D state and lost to the cooling process. 
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To measure the hyperfirie splitting of either 5D state, the ion can be initially placed in 
the (F g = 2, Mp = 0) Zeeman sublevel of the 651/2 ground state by optical pumping with 
a polarized 493 nm beam. As shown in Fig. [2J in the case of the 5D 5 / 2 measurement the 
ion is then transferred to a particular hyperfine sublevel (F, Mp) of 5D 5 / 2 by applying a 
pulse of resonant 1762 nm light. RF field coils are then turned on for driving a AF = ±1 
transition, after which a second pulse of the 1762 nm laser will transfer from (F, Mp) back 
to the ground state. The ion therefore ends up in the ground state if there was no hyperfine 
transition, and in 5D 5 / 2 if the RF is on resonance and the hyperfine transition was successful. 
In the former case there will be fluorescence when the ion is illuminated by the 493 nm/650 
nm lasers while in the latter case there will not be fluorescence; the ion is 'shelved' in the 
5D 5 / 2 state. The process is repeated for a range of RF frequencies and a hyperfine transition 
resonance curve is acquired. For the 5Ds/ 2 measurement, the same procedure is followed 
with 2051 nm resonant light to populate 5D 3 / 2 sublevels, but an extra step is needed at the 
end - the shelving of the 65*1/2 state population to the 5D 5 / 2 state. 

TABLE III: Zeeman splitting coefficients defined in Eq. (|18p . 
Mp 01 (MHz/G) a 2 (MHz/G) a 3 (MHz 2 /G 2 ) 



3 


3.359 


0.2099 


0.7052 


2 


2.239 


0.1400 


1.153 


1 


1.120 


0.06998 


1.422 











1.511 


-1 


-1.120 


-0.06998 


1.422 


-2 


-2.239 


-0.1400 


1.153 


-3 


-3.359 


-0.2099 


0.7052 



In the previous measurements of Zeeman transitions among 5D^/ 2 sublevels, sensitivities 
of a few Hz were achieved [ijj] , limited by incompletely shielded magnetic field fluctuations. 
In the hyperfine measurements proposed here, this source of broadening can be eliminated by 
using transitions that have a weak magnetic field dependence, such as Mp = — > M F = 0. 
In practice, it is often easiest to perform laser cooling in nonzero magnetic fields (typically 
B = 1G) to avoid formation of inefficiently cooled dark states. Such modest magnetic 
fields will introduce only a small B 2 dependence in the 0^0 transition when the energy 
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separations between hyperfine levels are sufficiently large. In fact the Zeeman effect for 
B = 1G can indeed be considered a small perturbation to the hyperfine splitting for all the 
5D 3 /2 states and for the F = 1 and F = 2 states of the 5D 5 / 2 manifold. 

However, for the 5D 5 / 2 F = 3 and F = 4 states the situation is more complicated 
because the hyperfine splitting between these states is small, = 0.49 MHz [26], so the Zeeman 
Hamiltonian cannot be considered a weak perturbation to the hyperfine Hamiltonian. The 
Zeeman and hyperfine Hamiltonians must therefore be treated on an equal footing. The 
manifold of F = 3 and F = 4 states are then all degenerate in zeroth order, and the 
first-order energy shifts due to Hzeeman + Hh yP are obtained by diagonalizing the matrix of 
H zeeman + Hh yp within the manifold of F = 3 and F = 4. In the \IJFMp) basis, the matrix 
is 2x2 block diagonal, with each block having a given value of Mp. The eigenvalues are then 
given by 

E± = \Ez + ai B±^ ±E% + a 2 E 3 B + a 3 B* (17) 

where E + reduces to the F = 3 hyperfine energy E 3 for B = and E_ reduces to the F = 4 
hyperfine energy, which is here taken as zero. The values of ai, a 2 , and 03 are given in Table 
III. These are related to the matrix elements of the Zeeman Hamiltonian H' as follows 

ai B = - ( (4, M F \H'\A, Mp) + (3, M F \H'\3, Mp) ) , 
a 2 B = - l -{ (4, M F \H'\4:, Mp) - (3, M F \H'\3, Mp) ) , 
a 3 B 2 = - ( (4, M F \H'\4, Mp) - (3, Af F |#'|3, Mp) f + ( (4, M F \H'\3, Mp) f . (18) 

A portion of the graph of the energy levels as a function of B is given in Fig. [3l Note 
that for even small magnetic fields there is a great deal of mixing between the F = 3 and 
F = 4 states, and there are two energy levels, originating from F = 3, Mp = — 1 and 
F — 4, Mp = 1, that are almost field-independent for B in the range of ~ 0.6G - 1.2G. 
Each of these levels can be connected to the F = 2, Mp = level by an RF transition, 
with very weak dependence on B. Likewise, the desired zero-field hyperfine intervals can be 
extracted from these RF measurements using only a relatively low resolution determination 
of B by a field-dependent Zeeman resonance. Thus measurement of all hyperfine intervals 
to 0.1 Hz seems feasible. 

We have shown that a simultaneous measurement of the hyperfine splittings in the 5D 3 / 2 
and the 5D 5 / 2 fine structure levels of Ba + allows one to unambiguously extract the value 
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FIG. 3: D 5 / 2 F = 3 and F = 4 hyperfine Zeeman levels. Energies are measured relative to the 
F = 4 energy at zero magnetic field. Very small mixing of the F = 3 state with the F = 2 state 
has been neglected. 

for the nuclear magnetic octupole moment. We performed the ah initio calculations of the 
relevant matrix elements in the framework of relativistic many-body perturbation theory, 
analyzing the first- and the second-order corrections to the hyperfine energy levels. We 
have also outlined an experimental procedure for measuring the hyperfine intervals to the 
required accuracy with single trapped Ba + ions. 
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